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Electron states are studied for quantum dots in a strained Si quantum well, taking into account 
both valley and orbital physics. Realistic geometries are considered, including circular and elliptical 
dot shapes, parallel and perpendicular magnetic fields, and (most importantly for valley coupling) 
the small local tilt of the quantum well interface away from the crystallographic axes. In absence 
of a tilt, valley splitting occurs only between pairs of states with the same orbital quantum num- 
bers. However, tilting is ubiquitous in conventional silicon heterostructures, leading to valley-orbit 
coupling. In this context, "valley splitting" is no longer a well defined concept, and the quantity of 
merit for qubit applications becomes the ground state gap. For typical dots used as qubits, a rich 
energy spectrum emerges, as a function of magnetic field, tilt angle, and orbital quantum number. 
Numerical and analytical solutions are obtained for the ground state gap and for the mixing fraction 
between the ground and excited states. This mixing can lead to valley scattering, decoherence, and 
leakage for Si spin qubits. 
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I. INTRODUCTION 

In the past decade, much progress has been made in 
the development of quantum dot devices for quantum 
information processing, spurred by the spin qubit pro- 
posal of Loss and DiVincenzo.'i' In GaAs devices, recent 
achievements have included single electron occupation of 
a quantum dot, Impulsed gating and single shot readout, SI 
and qubit rotations about an arbitrary axis.l^ However, 
it was noted early on that materials properties in silicon 
are potentially well suited for quantum computing.^ Pre- 
dominant among these are weak spin-orbit coupling and 
the ready availability of a nuclear spin-0 isotope. 

Several proposals were pu t forth for silicon-based 
quantum dot spin qubits.'^'I^ Many of the technologi- 
cal achievements made in GaAs have now been realized 
in Si. Recent experimental advances include controlled 
tunnel barriers in modulation dopedj^^^lSI degener ately 
doped,!^ and undoped enhancement-mode devices, '^^'^^ 
coupli ng betw een a quantum dot and a nearby charge 
sensor jHHSIiSl control led t unnel coupling between two 
sides of a double dot ]^°* ^^i single electron occupation of 
a quantum dotj^-- and observations of coherent quantum 
phenomena involving spins .12211211 

On the other hand, it is known that the multi-valley 
structure of the conduction band in Si may pose a chal- 
lenge for quantum computing.l23 In the context of quan- 
tum dot spin qubits in strained silicon heterostructures, 
it is necessary for the two-fold degeneracy of the two- 
dimensional electron gas to be lifted sufficiently that 
spins, not valleys, form the relevant two-level system.l^Sl 

The magnitude of the valley splitting has now been 
measured in strained Si quantum wells at both high 
fields^^ and low fields. '^"'^ In certain cases, the mea- 
sured values are 1-2 orders of magnitude smaller than ex- 
pected from theoretical calculations. '^^"^ Such splitting 
is too small for robust quantum computing applications. 
It appears that the unexpected suppression is caused by 
atomic-scale steps at the quantum well interface, as orig- 



inally suggested by Ando.l^ In this picture, the steps 
can arise from the underlying miscut in the growth sub- 
strate, which is transferred to the interface by conformal 
epitaxial deposition. Alternatively, they can be a conse- 
quence of strain-induced roughness. In either case, when 
an electronic wavefunction extends laterally over many 
steps, the valley splitting cannot be maximized on ev- 
ery step, leading to an overall suppression of the valley 
splitting. Since the valleys states include a phase factor 
whose value depends on the position of the interface,!^ 
the suppression of valley splitting can be interm-eted as 
an interference effect caused by multiple steps.l^ 

It should therefore be possible to overcome the sup- 
pression of valley splitting by laterally confining the 
wavefunction to a small number of steps. This proposal 
was put to the test in recent experiments by probing 
the valley splitting under strong electrostatic or mag netic 
confinement, in a quantum point contact geometry-l^SSD 
Under such conditions, the valley splitting was restored 
to its expected theoretical value. By the same token, a 
smooth interface could partially explain the exception- 
ally large valley splitting observed at a Si/Si02 inversion 
layer .ESl 

For quantum computing applications, we are most in- 
terested in quantum dot devices, where the confinement 
is controlled by the top-gate geometry and the applica- 
tion of gate voltages and magnetic fields. Recent experi- 
ments have shown that excitation energies in few-electron 
quantum dots may be of order 0.1-0.3 meV, which pro- 
vides an estimate for the valley splitting.^'* From a quan- 
tum computing perspective, the size of this splitting is 
encouraging, since it is larger than both the typical elec- 
tron temperature in a dilution refrigerator 150 mK), 
and the Zeeman spin splitting in fields up to about 2 T. 
(Note that valley splitting and Zeeman splitting are both 
functions of the magnetic field.) The magnitude of the 
valley splitting therefore appears sufficient for quantum 
computing in many cases. However, it is essential to un- 
derstand the complicated dependence of valley splitting 
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on the shape of the quantum dot and on the materials 
parameters associated with siheon heterostructures. In 
the context of decoherence, it is important to understand 
how vahey-orbit couphng causes quantum dot orbitals to 
hybridize. The latter effect plays a role in valley scatter- 
ing, and can allow the spin qubit to leak into a higher 
dimensional spin- valley Hilbert space. 

There are many possible approaches to modeling sil- 
icon quantum dots, ranging from highly simplified ge- 
ometries, such as circular quantum dots on a flat quan- 
tum well interface, to realistic three-dimensional struc- 
tures with top-gates and disordered interfaces. While 
the simple geometries may be treated analytically, the 
more realistic geometries require numerical methods. In 
this paper, we take a middle road by solving a somewhat 
complicated model that still allows analytical solutions. 
Specifically, we study inter-orbital valley coupling in a 
quantum dot with an elliptical shape and arbitrary step 
orientation, in parallel or perpendicular magnetic fields. 
In a previous numerical analysis, we also included the 
effects of interfacial disorder in circular dots,'^ although 
we did not account for inter-orbital valley coupling. In 
that case, we observed that interfacial step disorder has 
a strong moderating effect on the suppression of valley 
splitting. Because disorder is not easy to accommodate 
in our analytical treatment, we focus here on a uniform 
interface. However, we do allow for an effective tilt angle, 
which approximates some aspects of disorder. 

The paper is outlined as follows. In Sec. |ll] we de- 
scribe our theoretical method and procedure. In Sec. |III[ 
we consider a quantum dot geometry with the magnetic 
field oriented parallel to the sample surface. General so- 
lutions are obtained for the valley coupling in an ellipti- 
cal quantum dot, whose orientation with respect to the 
steps is arbitrary. We solve several important limiting 
cases, including high aspect ratios and circular quantum 
dots. In Sec. |IV[ we study circular and elliptical quan- 
tum dots in a perpendicular field geometry. In Sec. |Vj 
we summarize our main results and provide an intuitive 
explanation for the origin of valley-orbit coupling. We 
discuss some consequences of valley-orbit coupling and 
why valley splitting is not an entirely well-defined quan- 
tity. We also discuss materials parameters that affect 
the magnitude of the valley splitting, such as disorder 
at the interface, and we suggest future experiments to 
characterize these materials properties. In the appendix, 
we provide a discussion of the important approximations 
used in this paper, including the treatment of the the 
interfacial tilt as smooth and uniform. 



II. VALLEY-ORBIT COUPLING 
A. Effective mass formalism 

We build on the work of Ref. 36 to develop an effec- 
tive mass theory of valley coupling in a quantum dot. 
In contrast with most previous theories, we now include 



inter-orbital couplings. These have a significant effect on 
the valley splitting in typical qubit devices, and they are 
essential for understanding valley scattering. 

In this theory, we treat the valley coupling as a per- 
turbation. Thus, at zeroth order, valleys do not enter 
the analysis except through the anisotropic effective mass 
tensor. The resulting equations have a simple form for 
inversion layer^^and strained silicon heterostructures,'^ 
since only z valleys are low enough in energy to play a 
role in the calculations. In this case, there is only one 
envelope function equation: the effective zeroth order 
Hamiltonian, 
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-fVQw(2') + ^QD(x',2/') 



(1) 



In Eq. ([I]), we refer to two different coordinate systems. 
The unprimed coordinate system (x, y, z) is aligned with 
the crystallographic axes, while the primed coordinate 
system (x', y', z') is aligned with the growth axis, where 
z' is perpendicular to the plane of the quantum well. For 
the z valleys, the transverse effective mass is given by 
nix = TTiy = rrit 0.19tooi while the longitudinal effec- 
tive mass is given by niz = mi 0.92too. The confine- 
ment potentials include Vqw, the vertical quantum well 
potential, and Vqd, the lateral quantum dot potential. 
The eigenstates of Hq are denoted as envelope functions, 
F„(r'), with the discrete orbital index n. Since Hq does 
not depend on the individual z valleys, the energy eigen- 
values £„ are doubly-degenerate. 

At first order in the perturbation theory, we introduce 
a valleycoupling potential T4 . Following Fritzschc'^ and 
Twose,'^we can obtain a set of coupled equations for the 
two- valley system: 

J2 E «n,„,.e""=«-[i7o + K(r') - E^]F„ir') = 0. (2) 

n v—dzl 

Here, the index v ~ ±1 refers to the valley centered at 
ky = vkQZ, where /cq — 0.82(27r/a) and a = 5.43 A is the 
Si cubic unit cell dimension. The interaction Vy mixes 
unperturbed orbitals with different quantum numbers n, 
so Eq. ([2| must include a sum over the index n. The re- 
sulting eigenstates of the perturbation Hamiltonian are 
indexed by a new label, n, with corresponding eigenval- 
ues En and eigenvectors an, whose components span the 
unperturbed basis set {|n) (g) \v)}. Due to valley-orbit 
mixing, it is not appropriate to distinguish between or- 
bital and valley quantum numbers in realistic quantum 
dots, except in certain limiting cases. The index n is 
therefore a combined valley-orbital label. 

We note that Umklapp terms are absent in the sum- 
mation of Eq. ([2|. In the context of valley coupling 
by a shallow donor, the omission of Umklapp processes 
is known to give numerical errors in the calculation 
of binding energies .'^^'^ While these might be resolved 
by going beyond conventional effective mass theory,!^ 
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the conventional theory does provide acceptable results 
and intuitive explanations that are in agreement with 
more detailed analyses, both for single-valley and multi- 
valley semiconductors!^ For quantum wells, the two- 
valley treatment of Eq. ([2| provides remarkable agree- 
ment with more sophisticated theories P^I^fil This is not 
surprising - one would expect the Umklapp errors to be 
suppressed in quantum wells, as compared to donors, be- 
cause the fourier transformed potential Vqwik) decays 
more quickly in reciprocal space. Based on the overall 
agreement between effective mass and more detailed the- 
ories, we will not pursue Umklapp processes here. 

Each term in Eq. ^ corresponds to a different unper- 
turbed valley orbital, identified by indices v and n: 



(r>,i;) =e''='"^^F„(r'). 



(3) 



It is also convenient to define a special index referring to 
the opposite valley, v = —v, such that 

(r>,i;) = e-"=«'^^F„(r'). (4) 



As discussed elsewherej ^^ * '^^ * ^^ ^ the valley coupling po- 
tential takes the form of a (^-function in the effective mass 
theory: 



K(r') = Vv6{z' - zo), 



(5) 



where z' — zq is the position of the sharp quantum well 
interface. In principle, one (5-function potential should 
be included for each interface of the quantum well. How- 
ever, in a typical modulation-doped device, the envelope 
function is centered near the top interface (closest to the 
doping layer) , and is exponentially suppressed at the bot- 
tom interface. Therefore, we include only one (5-function, 
for the top interface. (If necessary, a second interface 
could easily be included in the theory.) 

Equation ([S]) describes a geometry where the interface 
occurs on a horizontal plane in the (x',?/',z') coordinate 
system. This assumes that the stepped interface can be 
approximated as uniformly tilted. The validity of such 
an approximation is discussed in detail in the Appendix. 
In brief, since a single-atom step of height a/4 = 1.4 A is 
much smaller than any effective mass length scale, the 
continuum approximation is well justified in most situa- 
tions. Additionally, the assumption of a smooth interface 
is well-justified, locally, in the absence of significant step 
bunching. 

When we compute matrix elements, as described be- 
low, it is more convenient to transform to the crystallo- 
graphic (cc, y, z) coordinate system: 



K(r) = Vv5[z - Zi{x,y)), 



(6) 



where 



Zi{x,y) = zq — X tan I? cos — y tan 1? sin iy9 

~ zq — xd cos ip — yd sin ip. (7) 

Here, d is the interfacial tilt angle, which is typically 
small. For example, samples grown on commercial sub- 
strates may have 1) < 2° , corresponding to a step width of 



3.9 nm. ip is the tilt direction in the x-y plane, with ip = 
corresponding to a downhill slope in the x direction. For 
a quantum dot covering many steps, the valley coupling 
should depend very weakly on the relative position, zq, 
as we have verified numerically. In the opposite limit of 
very few steps, we should treat the steps explicitly, as de- 
scribed in the appendix. For the present analysis, where 
we treat the interface as smoothly tilted, it is appropriate 
to choose the value of zq to simplify our calculations. 

The parameter should be viewed as an effective tilt 
angle. (See appendix.) Normal growth dynamics tend 
to produce interfaces that are disordered,'33 so the local 
tilt angle can differ from the global tilt. As discussed be- 
low, lower ground state energies can be attained at flat 
interfaces. Confined electrons will therefore adjust their 
center positions slightly, to seek out locally flat regions. 
A more careful treatment of this phenomenon, and disor- 
der in general, lies outside the scope of the present work. 
However, we can capture some of the relevant physics 
simply by treating as an effective, locally-averaged tilt 
angle. 

To solve Eq. |2]), we proceed as in Ref. |3S1 we left- 
multiply by e~™ '^"'^ F^,{r'), integrate over all space, and 
remove any resulting small terms. In particular, we drop 
terms in the integrand containing fast oscillations.^^ In 
this way, we obtain the following expression for the ma- 
trix elements of the perturbation Hamiltonian: 

{n\v'\H\n,v) (8) 
= [ e*(^-^')'=°^F,:,(r')i?F„(r')dr', 



where we define H = Hq + Vy. 

After solving for the valley orbitals, we may finally 
combine the atomic and confinement scale components 
of the wavefunction, following the prescription given by 
KohnP 

*a(r) ~ ^[a;i,„,_ie-"=«^u_fe„(r) (9) 

n 

+a^i,„,+ie+''=°^M+fc„(r)]F„(r'). 

Here, e^'^'^''^u±kg{v) are the bulk Bloch functions for the 
two valley minima. Note that the periodic functions 
u±ko (r) do not play a role in our calculations of the valley 
orbitals. The atomic scale details of the valley coupling 
are captured in the coupling parameter Vy. 



B. Selection rules and broken symmetries 

We evaluate the matrix elements of Eq. ([8| as follows. 
For the case v' = u, we find 



(n', vlHln, v) 



(10) 



£n-^Vv / |F„(x',?/',zo)| dx'dy 
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where we have used Eq. ([s]). For every quantum dot 
geometry studied in this paper, we observe the following 
form of separability for the envelope function: 



|i^„(a;',y',z')|^|/„(a;',y')IC(^')- 



(11) 



Eq. ( 10 ) then reduces to 

(n', v\H\n, v) = Sn,n' [e™ + Vv('^{zo)] , (12) 

when the envelope functions are properly normalized. 
The case v' — v can be expressed as 

{n',v\H\n,v}=v,e{zo)e^'^''"'" (13) 

--2ivkQ [x 1? cos ip-^y^'d sin if] 

xf:,{x\y')U{x',y')dx'dy'. 

Here, we have made use of Eq. ([t]), and the fact that 
{x, y) ~ (x', y') when ■& ~0. Note that we have dropped 
the terms {n' ,v\H(i\n,v) because they contain fast os- 
cillations and are therefore small. Since all coordinate 



variables in Eq. ( 13 1 are primed, we henceforth drop the 



prime notation for simplicity, except where noted. 

When v' — V and = 0, Eq. (13) simplifies, due to 
the orthonormality of the valley orbitals. The result is 
proportional to However, when i9 7^ 0, the phase 

factors inside the integral are non-vanishing, and inter- 
orbital couplings arise. In this way, we see that the tilt 
of the quantum well is directly responsible for the off- 
diagonal valley-orbit coupling. Since roughness is ubiq- 
uitous in realistic devices, we must always expect some 
degree of inter-orbital valley coupling. 

The emergence of inter-orbital coupling has further 
interesting consequences. First, it causes broken time- 
reversal symmetry. When = 0, the perturbation pro- 
duces only intra-orbital terms; the valley coupling is iden- 
tical for every orbital, and its phase factor is trivial, since 
it can be eliminated by setting zq = 0. But when 7^ 0, 
the inter-orbital couplings are complex, and they can- 
not be trivially removed. The perturbation Hamiltonian 
is intrinsically complex, and time reversal symmetry is 
broken.^" Similar behavior is observed in the multi-valley 
graphene system, with important implications for weak 
localization.l^i^l^ Second, the valley index can no longer 
be treated as a good quantum number. When = 0, the 
valley labels are universal, in the sense that they have 
the same meaning for every orbital. In a square well, for 
example, the wavefunctions can be unambiguously char- 
acterized as odd or even.^^ However, when ^ 0, there 
is no universal labeling scheme. Valley-orbit coupling 
causes level mixing, which differs from orbital to orbital. 



C. Summary of approximations and their 
consequences 

The effective mass approximation forms the basis of 
our analysis, as it allows us to analyze quantum dots in 



terms of their envelope functions rather than their wave- 
functions. The approximation is a natural one, based on 
the large separation of lengths scales!^ atomic (^0.1 nm) 
vs. confinement (10-100 nm). One could expect the ap- 
proximation to break down when the confinement poten- 
tial is sharp, since the two length scales are then com- 
parable. However, the effective mass approximation re- 
mains viable for many materials, including silicon, even 
near an abrupt interface.!^ This is because the valleys 
are roughly parabolic in the regions of the Brillouin zone 
where the wavefunction is concentrated. 

Valley coupling is therefore a higher order effect in the 
effective mass approximation. A confinement potential 
with components at large k will couple the valleys in sili- 
con, leading to an envelope function equation of the form 
of Eq. ([2]), where the leading order approximation for the 
coupling takes the form of a 5- function, as in Eq. (|6|. 
When the interface is aligned with the crystallographic 
axes, the interface position Zi is a constant. This leads 
to our first selection rule, that a sharp interface always 
causes inter-valley coupling within a given orbital. We 
note however, that even if the effective mass approxima- 
tion were to break down near a sharp interface, the valley 
degeneracy would still be lifted, as confirmed by atom- 
istic theories.lS^l Comparison between atomistic and effec- 
tive mass theories generally confirms the effective mass 
approach. ^''^^ In the absence of any interfacial tilt, the 
valley coupling is strictly intra-orbital. This is a con- 
sequence of the assumed separation of variables in the 
confinement potential, which forms an excellent approx- 
imation in quasi-2D devices.'^ 

When the interface is not aligned with the crystal- 
lographic axes, we have chosen to treat the interface 
as smoothly tilted, since single atomic steps are much 
smaller than effective mass length scales. The approx- 
imation leads immediately to our second selection rule, 
that a tilted interface causes inter-orbital coupling be- 
tween opposing valleys, v and v. Of course, the smooth 
interface approximation will not be valid if the misalign- 
ment is highly non-uniform, as in the case of multiple 
step bunching. In the appendix, we consider several sce- 
narios where the smooth interface approximation breaks 
down, most notably in cases involving wide steps. For 
very wide steps, the smooth approximation leads to a 
very weak valley coupling. A more accurate treatment, 
taking into account the relaxation mechanism described 
in appendix A.l, would suppress the valley coupling even 
beyond the predictions of the present theory. 



III. PARALLEL FIELD GEOMETRY 

We now analyze some specific geometries. We first 
consider the case where the magnetic field is oriented 
parallel to the sample. Because the x and y directions 
are now inequivalent, even in the absence of interfacial 
steps, the angular momentum quantum number I is not 
a good quantum number. 
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We initially treat the problem as an anisotropic 3D 
simple harmonic oscillator, with lateral confinement fre- 
quencies LJx and LOy. In the most general case, we as- 
sume the oscillator frequencies are unequal: ujx 7^ '^y 
The quantum well confinement potential is also taken to 
be parabolic, with lOz 3> LOx,^y Because of the strong 
confinement along z, the quantum dot wavefunction is 
approximately separable in the variable z, despite the 
presence of a magnetic field. However, we will find that 
the wavefunction acquires an intrinsic phase when the 
interface is tilted, due to valley coupling. Since the mag- 
netic field also produces a non-trivial phase, which could 
potentially lead to interference effects, it is important to 
work through the details of the problem without assum- 
ing separability, a priori. 

We take the magnetic field to be oriented along one 
of the major elliptical axes (y), otherwise the problem 
becomes intractable. Adopting the symmetric gauge 
A — [z, 0, —x)B/2 for the vector potential, the envelope 
Hamiltonian becomes 



^0= E 



1 

2™, 



d 



2 99 



(14) 

The eigenvalue problem is separable in the variable y, 
leading to solutions of the form 



F„^,ny{r) = fnAx,z)g„ (y), 



(15) 



where and riy are non-negative integers. The solutions 
for 5„ are given by 



,„„(,) = (2-n,!)-/^(";|^)^^%-(- 



mtujy 



(16) 



where Hn{x) is a Hermite polynomiaLl^ 

The magnetic field couples the x and z variables in 
Eq. ( |14[ ), making the problem somewhat more compli- 
cated. However, by rescaling the coordinate variables as 



X = [mi/mtY^^x, 
z = [mtlmiYl^z, 



(17) 
(18) 



the x-z Hamiltonian is brought into a more tractable 
form: 



1 

2^ 



-ihV + eA 



/ 2-2 I 2~2\ 



(19) 



where /i = y/rritmi, and we define a scaled vector poten- 
tial, A = (z,0,-x)B/2. 

In principle, the uJz oscillator possesses a distinct quan- 
tum number, Uz- However, we restrict our investiga- 
tion to the first subband approximation (n^ — 0), since 
UJz 3> iOx- We also assume that the quantum well 
confinement is stronger than the magnetic confinement, 
LUz S> uJc, as consistent with typical laboratory condi- 
tions. (Here, ojc = e\B\/mt is the cyclotron frequency.) 
Although valley coupling can cause mixing of the higher 
subbands, the relevant energy scales are very large, and 
the couplings are small. We ignore such mixing here. 



An exact solution for Eq. ( 19 ) was first obtained in 



Ref. [571 For the problem considered here, the results can 
be expressed as 



fn^{x,z) = Nn^ exp 



iJxz MiijJi{mi/mty/^z^ + 2iLAIiUJiM2CjJ2XZ + M2UJ2{mt/miy/^x^ 



2h{L'^MiUJiM2UJ2 + 1) 
' ^2M20J2/n{LMiUJi{mi/mt)^''^z - i{mt / mif / ^ x) 

L'^MiUJiM2L02 + 1 



(20) 



where 



(21) 



Ml, 



(23) 



2^1^ (ujI + gjg + a.g)2 _ 4^2^2 



(22) 



1 

7i 



2 I 2 I 2 



(24) 



±v/(w2+cj2+w2)2_4^2^2 



,1/2 



6 



Cl 



y/MiUJiM2UJ2 

n'^h{n^y.{L^MiUJiM2UJ2 + 1) 



r (4G)'/2r(^) (Zeven), 



I 



(/ odd), 



and 



L^MiUJiM2UJ2 + 1 



- 1. 



1/2 

, (25) 



(26) 



(27) 



Here, r(p) is the Gamma function!^ The energy eigen- 
values corresponding to Fn^^uy are 



We now apply the limit lo^ S> uJj.,uJc, obtaining 

,2 



U2 ^ UJx+ 0[uj^^]. 



(29) 
(30) 



Including the leading order correction in (u)c/^z), the 
quantum dot energy is then given by 



2 r^^^ 



Tlx + M hhJj, + (uy + ^ 1 fiw,, 



(31) 



with the corresponding x-z eigenfunction 

mtUJx 2 , i^/mtmiUJc miuj^ 2 
X exp <, H — - — — xz — — z 



(32) 



2h 



2h 



2h 



mtUJx 



Here, we have made use of a Hermite polynomial identity: 

2n/2„ ( ^_+y\ A n\ 



E 



V2 ; ^^k\{n~k)\ 



Hk{x)Hn-k{y)- (33) 



We find that the main effect of the magnetic field is 
to introduce a phase factor into the otherwise separable 
anisotropic oscillator solutions. This phase plays no role 
when it is included in valley coupling calculations, since 
it is independent of n^. The conventional wisdom that 
quantum dot wavefunctions are separable in the paral- 
lel field geometry is therefore also valid in the context of 
valley coupling. 




FIG. 1: (Color online) Energy levels for an elliptical quan- 
tum dot in the parallel field geometry, including valley cou- 
pling. Dot parameters are appropriate for a single-electron 
spin qubit: hujx ~ 1.5 meV, hujy — 0.75 meV, w^Co = 0.5 meV. 
Uniform, lateral step orientation is set to = 7r/4. Blue 
dot-dashed lines show the analytical approximation, obtained 
by r estri cting the analysis to the lowest two orbital levels in 
Eq. ( |46[ ). The subband energy shift (a constant) has been 
dropped here, for simplicity. 



A. Flat Quantum Well 



We 



now 



iKn'yV'lHln^ny. 



compute the matrix elements 
yv) for a flat quantum well, treating 
the valley coupling potential Vy as a, perturbation. 
When v' = V, only diagonal matrix elements survive, 
as explained in Sec. |nB[ Note that we have chosen the 
center of the harmonic oscillator to be at the origin. 
Defining the top of the quantum well as z — zq, we 
obtain 



{nxnyv\H\nxnyv) = e,. 



with 



Co(zo) = (m^c.,M)l/4e-("'"^/2'')-o^ 



(34) 



(35) 



This final expression was obtained by treating the vertical 
confinement potential as harmonic. 



The approximate separability of variables in Eq. ( 32 ) 
suggests that hcuz and Co could be replaced by more 
physically motivated quantities, such as eigenstates of 
a finite triangular well. Indeed, throughout this work, 
we will view VvCo ^ scaling factor for the valley cou- 
pling, whose value can be measured experimentally and 
inserted, phenomenologically, into the theoretical expres- 
sions. When microscopic disorder is present at the quan- 
tum well interface and significant details about the val- 
ley coupling are not known precisely, such an approach 
becomes quite practical. Consistent with this view, we 
henceforth drop the argument zq in Cq. 
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elements 
J - the eigen- 

functions restricts the off-diagonal couplings to cases 



For the ofF-diagonal matrix 
{n'^n' v\H\nx'nyv) , the orthogonality of 



where 
(Recall that 'd' 
result 



and 



as discussed in Sec. 



II B 



in this section.) This gives the trivia 



2 2ifco zov 



(36) 



The two-fold degeneracy of each orbital level is lifted by 
a fixed amount, independent of the quantum state. The 
ground state gap for each pair of states is given by Eg — 
2vyC,Q. This is the conventional "valley splitting." 



B. Tilted Quantum Well 

We turn to the more realistic case of a quantum well 
tilted downwards at an angle i?, towards the direction if 



in the x-y plane, as described in Eq. Q. When v' = v, as 
in Eq. ( 10 ), then to leading order in the small parameter 



{ujcl^z)i the Hamiltonian matrix elements are real and 
exclusively diagonal. They are identical to their counter- 
parts for a flat quantum well. However when v' — v, the 
valley interaction induces nontrivial couplings between 
the quantum dot levels, leading to new avoided cross- 
ings in the energy spectrum, and to fundamentally new 
physics. 

The coupling matrix elements are solved using an ex- 
plicit expansion for the Hermite polynomials. 



L"/2J , . 



(37) 



where the floor function [mj is defined as the greatest 
integer less than or equal to m, and n\ is the factorial 
function with 0! = 1. We obtain 



7r2;i22».+<+»«+»;(„^)!(„g!(„^)!(„g! 



1/2 



dxdy exp 



mtUJx 2 mtujy 2 

— fc — ^ — i~y 



exp [—2ik{)v{x'd cos tp -\- yd sin </?)] 



mtUJx 



mtUJy 



rutuiy 



= t;„Co'e'^'°^''VK)!K)!K)!K)!exp 

Ln,/2J L";/2J L"„/2J L»;/2J 



hkQ'd'^ f cos^ ip sin ip 



mt 



(38) 



— 
iV2j 



X E E E E [^!(O!(«.-20!«-2O!j!(/)!(n,-2j)!(n;-2/)!]-' 

/=0 l'=0 j=0 j'=0 



X H„ 



,+<-2(/+i') 



mtUJx 



-kQ-d cos if 



v+n'y-2{j+j') 



rritujy 



-k^d sin Lp 



This result is separable with respect to the x and y in- 
dices, {ri3;,n^} and {rty,nj^}, indicating that mixing oc- 
curs among x orbitals or y orbitals, with no x-y cross- 
coupling. At this point, we invoke our freedom to set 
zq = 0, as discussed in Sec. |II A[ 



The couplings appearing in Eq. (38) suggest a rather 



complicated dependence of the valley splitting on the tilt 
angle The energy spectrum for typical dot parameters 
is shown in Fig. [l] The different asymptotic behaviors 
can be understood as follows. In the limit — > 0, rota- 



tional symmetry is recovered (see Sec. II B), and valley 
coupling does not mix the orbital levels. In this case, the 
valley states in each orbital level are split by the same 
amount. However, for realistic quantum dot diameters 
(~ 50— 100 nm), the orbital spacing can be smaller than 



the single-particle valley splitting (^ 0.3 — 0.5 meV), as 
in Fig. [l] In this case, pairs of states with the same or- 
bital quantum numbers are difficult to identify, due to 
the presence of valleys. 

When I? > 0, valley-orbit coupling mixes the unper- 
turbed orbitals, leading to avoided crossings. The degree 
of mixing affects the valley scattering. We can quantify 
this effect by defining the mixing fraction J^a^o, which 
describes the projection of the perturbed ground state 
{ii = 0) onto the unperturbed excited states (n > 0): 



E 

n>0,D=±l 



ah. 



(39) 



Some typical results for the mixing fraction are shown 
in Fig. [2j We see that J^h^o vanishes when i9 — >• 0, as 
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FIG. 2: (Color online) Mixing fraction, Eq. ( [39| ), due to inter- 
orbital valley coupling. Results are shown for the same quan- 
tum dot parameters used in Fig. [T] The solid black line is 
the projection of the perturbed ground state (n — 0) onto 
the unperturbed excited states (n > 0). The red dashed line 
shows the main contribution to J-n^o coming from just the 
first excited orbital {n — 1). Blue dot-dashed line shows the 



corresponding analytical result of Eq. ( 46 1 , as obtained within 
a restricted subspace. Inset: the projection of the first excited 
state (n — 1) onto the unperturbed excited states (n > 0). 



explained above, and goes through a maximum value. In 
the limit ^ » 1, the valley coupling is suppressed for 
all orbital states. (See Fig. [T]) This is a consequence of 
destructi ve va lley interference, caused by multiple step 
coverage .l^iEII The mixing fraction is also suppressed for 
the same reason. In Sees. |III Cf E, we will obtain ana- 
lytical expressions for the valley-orbit coupling and the 
ground state gap, for certain quantum dot geometries. 

The first excited state experiences stronger valley-orbit 
coupling than the ground state, because of its proxim- 
ity to other nearby energy levels. This is demonstrated 
vividly in Fig. [T] In the two limits = and ^ 1, the 
orbital quantum numbers n^^y are good quantum num- 
bers. When I? = 0, the valley splitting is larger than 
the orbital spacing, so ny will typically vary on succes- 
sive energy levels (in this case, Uy = 0,1,0,..., from 
bottom to top). However, when ^ 1, the lowest pair 
of eigenstates becomes nearly degenerate, with the same 
quantum number, Uy = 0. Thus, the orbital character of 
the first excited state changes completely, from rij, = 1 to 
0, as a result of valley-orbit coupling. In the intermedi- 
ate region, i9 < 1, valley-orbit coupling causes a complete 
mixing of the unperturbed orbitals. The orbital crossover 
can also be observed in the inset of Fig. [2] which shows 
the mixing fraction for the first excited state. Clearly 
we cannot obtain exact expressions for the ground state 
gap or the mixing fraction. However, approximations are 
available in certain limiting cases, as described in the 
following sections. 



C. Limit: Small quantum dots 

This limit corresponds to the case where orbital split- 
tings are much larger than the valley coupling, leading to 
particularly simple approximations for the ground state 
gap. To leading order, we can treat the ground state by 
simply ignoring inter-orbital valley couplings. The quan- 
tum numbers = Uy = then remain good quantum 
numbers, even after applying the perturbation. 

We can express the Hamiltonian in the restricted man- 
ifold {|00+), |00— )}, where -I- and — refer to the un- 
perturbed valley quantum numbers. From Eq. (38) and 
Sec. |II B[ we obtain 



H : 

where we define 



£0 + Ao 
Ao?7 



Ao?7 
£0 + Ao 



£o 




Ao 


— ^vCqj 


V 


= exp 



TTlf 



COS^ (p 



sin^ (fi 



(40) 

(41) 
(42) 

(43) 



Note that the vertical confinement energy {h/2)(Ldz + 
bj^jlujz is a constant here; we set it to zero for simplicity. 

Diagonalization immediately gives the ground state 
gap, 



~ 2Ao77 
= 2t;„Coexp 



(44) 



rat 



cos^ Lp 



sin^ Lp 



The quantity 2Ao can be recognized as the theoretical 
maximum of the valley splitting, while 2A077 corresponds 
to the renormalized or suppressed valley splitting, which 
takes into account the steps at the interface. 

The previous approximation is most appropriate for 
ultra-small quantum dots in a 2DEG , with radii 
20 nm. Although such devices are challenging to fab- 
ricate at the present time, the main result of this section, 
Eq. (44), is more broadly applicable. This is particu- 



larly true when -d is so large that the valley coupling is 
strongly suppressed. We can use Eq. (44) to determine 
the crossover angle i?^, above which the valley coupling 
is exponentially suppressed: 



ds ^ \Jmtu:/hkl. 



(45) 



For a typical quantum dot energy of hu = 0.75 meV 
(see Fig. [I]) , this corresponds to a crossover angle Us — 
0.26°. Taking into account the simple harmonic oscillator 
radius, R = y''h/2mtuj, we obtain a relation between R 
and ^s, as shown in Fig. [3j 



D. Limit: High aspect ratio 

We now consider slightly larger quantum dots, of inter- 
est for quantum computing. In this regime, the orbital 
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Miscut tilt, i3 (degrees) 

FIG. 3: Crossover radius for the suppression of valley cou- 
pling caused by a tilted interface. For a quantum dot ly- 
ing below the curve, valley splitting is weakly suppressed or 
unsuppressed. Above the curve, valley splitting is strongly 
suppressed. 



and valley splittings are similar in magnitude. We further 
assume that one principle axis of the elliptical potential is 
much larger than the other (cOx ^ Wy). Inter-orbital cou- 
plings to the ground state are then dominated by a single 
excited state, Uy = 1. We find that this approximation is 
applicable over a fairly wide range. For example. Fig. [2] 
corresponds to a moderate aspect ratio of uj^/i-Oy ~ 2. 
The red dashed line shows results obtained by solving 
the Hamiltonian with many orbital states, but including 
just the contribution from the first excited state. The 
results arc satisfactory, indicating the accuracy of the 
approximation. 

For the analytical treatment described above, we con- 
sider the ordered manifold {|00-t-), |00-), |01-|-), |01-)}, 



(sa + Ao Ao?? il^nm 

Aqt] £o + Ao -iAoJyg 

iAom ei + Ao Aov{l - q^) 

-iAom Ao?7(l-g2) ei + Ac 



(46) 

Here, we use the previous definitions for eg, Ag, and rj, 
and we further define 



q = {koi3 sin (p) \J 2h/mtUJy . 



(47) 
(48) 



As before, the vertical confinement energy is a constant, 
and we set it to zero. 

Diagonalizing H gives the energy spectrum in this re- 
stricted manifold. Fig. [T] shows two two lowest energy 
levels obtained by this procedure (blue dot-dashed line). 
The results are satisfactory, even though the high aspect 
ratio limit is not well satisfied. In particular, the ground 
state gap is found to be quite accurate. 

One main advantage of considering restricted mani- 
folds is that they allow analytic results. Thus, the ground 
state gap can be computed, giving 



t Aam^ (49) 
\^[nLOy + Aar^{2-q^)Y + {2A^^)W 
\^[hLOy - Aor?(2 ~ g2)]2 + {2A^^)W- 



We can also compute the mixing fraction defined in 
Eq. pgl): 



2(Aor7)2g^ 

[hLOy + Ao77(2 - g2)]2 + (2Ao77)2g2 + [hwy + Aov{2 - q^)W[hwy + Aor/(2 - g2)]2 + (2Ao7?)2(j2 ' 



(50) 



In Fig. [2] we compare this analytical theory (dotted line) 
with the full numerical calculation. Again, we find that 
the restricted manifold provides a very reasonable ap- 
proximation. 



E. Limit: Circular dot 

Circular quantum dots are of interest for their sym- 
metry properties. We therefore introduce the angular 
momentum quantum numbers, n and I. In principle, a 
parallel magnetic field will break the x-y rotational sym- 
metry. However, for typical laboratory fields, we have 
shown that the only effect of a parallel field is to add a 



phase factor to the wavefunction. The phase factor does 
not appear in the Hamiltonian matrix elements. With 
this caveat, we can ignore the parallel magnetic field in 
the present analysis, and consider the wavefunction to be 
symmetric with respect to x and y. 

In Sec. |IV A[ below, we solve the circular quantum dot 
geometry in a perpendicular magnetic field. According to 
the previous discussion, we can obtain results appropriate 
for a parallel field from that analysis, by taking the limit 
-B — > 0. This amounts to making the substitution 



b/l% — > 2ujamt/h 



(51) 



in Eqs. (60l-(63). The corresponding energy levels are 
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given by 



where Aq is the same as before, and we now define 



huJz 



huJo{2n+\l\ + 1) 



(52) 



(parallel -B-field). 



Here, we have also included the leading contribution 
to the magnetic confinement energy, in analogy with 
Eq. ISTI). 



The energy spectrum of a circular quantum dot is 
qualitatively similar to an elliptical dot, and we do not 
show any plots here. However, some limiting behav- 
iors are of interest. In the small dot limit, we evaluate 
the perturbation Hamiltonian in the restricted manifold 
{\n,l,v)} = {|0,0,+),|0,0,-)}, similar to Sec. [!irC 
obtain the ground state gap 



We 



Eg ~ 2AoV, 



(53) 



rj = e ^ with q = k^d \/hJujQmt . (54) 



The result is equivalent to Eq. (|44|) in the isotropic limit 



ljJq, where we set tp — Q. 



We can also obtain analytical approximations for the 
circular quantum dot geometry when the orbital and 
valley splittings are similar in magnitude. We consider 
the manifold consisting of the ground state and the 
lowest excited orbitals. From Eq. (52 1, we note that 



±Z states are degenerate. The manifold of interest is 
therefore 6-dimensional, with the quantum numbers 
n — I — 0, ±1, and v = ±1. In the ordered basis 
{|0, 0, +), |0, 0, -), |0, 1, +), |0, 1, -), |0, -1, +), |0, -1, -)}, 
we obtain 



H 



Aor? £0,0 + Ao -iAoryg 

iAom £0,1 + Ao Ao?7(l-q2) 

-iAoi]q Ao?7(l-g2) eoi+Ao 

iAom -Aom^ £0,1 +Ao Aory(l-g2) 

V -iAom -Ao?7g2 o Ao77(l - q^) £o.i + Aq / 





-iAom 


-Aom^ 



iAom \ 


-Aow^ 




where we have used previous definitions, including Eq. (54 1 



Diagonalizing H gives the effective energy spectrum, and the resulting ground state gap, 
1 



Eg ~ 2Aor;g" 



^[hwo + 2Ao?7(l - q^)Y + 2(2Ao77g)2 - ^/[hujo - 2Aory(l - g2)]2 + 2{2Aomf 



We can also compute the mixing fraction: 



n=0 



{2Aom? 



[Hloo + 2Aor7(l - q^)] + '^{'^AomY + [^^0 + 2Ao77(l - 9^)] ^[huo + 2Ao7?(l - q^)f + 2(2Aow)' 

I 



(55) 



(56) 



(57) 



These results can be compared to an elliptical dot in the 
high aspect ratio limit (Sec 
ujo and if — )■ 



tutions ujy 



IIID[ ). Makin^the substi- 
7r/2 in Eqs 



consistent with our choice of lo^ 3> ojy 



(|49| and 
we obtain identi- 



cal results for Eg and J-n=o, up to leading order. Details 
of the energy spectrum begin to differ only for higher 
energy levels. 



IV. PERPENDICULAR FIELD GEOMETRY 

A. Circular dot 

In contrast with the previous section, the magnetic 
field plays a non-trivial role for valley coupling in the 
perpendicular field geometry. Here, we consider a single 
electron in a circularly-symmetric parabolic confinement 



potential, in a perpendicular magnetic field, B = Bi. 
(We continue to use unprinied coordinate notation.) The 
unperturbed Hamiltonian is given by 



^0= E 



1 

2m, 



n 2 



+Vqw(z) + ■7^'mtujl{x^ + y^), 



(58) 



where uo characterizes the parabolic potential. We adopt 
the symmetric gauge for the vector potential: A(r) = 
{—By,Bx,0). In this case, the eigenfunction function is 
separable: 



F{r) = f{x,y)C{z). 



(59) 



Appropriate solutions for the quantu m we ll envelope 
function, C(^); are discussed elsewhere !^^ * ^^ * As before. 
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we adopt the lowest subband approximation. We also 
choose the position of the interface, zq = 0, to simplify 
our calculations. 

The lateral eigen states f{x, y) are known as Fock- 
Darwin states.- The properly normalized Fock- 
Darwin functions can be expressed in terms of radial co- 
ordinates, {r cos 9 ^ r sin 6) = {x,y), as follows: 



fni{r,0) 



i\b 



27rZ|2l'l(7i+ |/|)! 



(60) 



xe e 



L\i\ 



2/2" 



Here, n and / are the orbital and angular momentum 
quantum numbers, respectively, n and I are both in- 
tegers, with n > 0. Additionally, Ib = \/f^/\^B\ is 
the quantum magnetic length scale, and L^^i^) ^ 
generalized Laguerre polynomial.!^ xhe quantity b — 



\/\ + {2ujq/u:cY is defined in terms of the cyclotron and 
confinement frequencies. As a consequence of the cir- 
cular symmetry, we can choose the step direction to lie 
along X without loss of generality, so that = 0. The 
Fock-Darwin energy levels are then given by 



em = ^ + ^[b{2n+\l\ + l) + l] 



(61) 



(perpendicular i?- field) , 

where we have included the subband energy for complete- 
ness, although it will be ignored in the following calcula- 
tions. 

Consistent with previous sections, the intra- valley per- 
turbation matrix elements are given by 

{n'J',v\H\n,l,v) = (5„^„'(S;,/' (£„/ -I- Ao), (62) 

while the inter-valley matrix elements are given by 



{n' ,1' ,v\Vy\n,l,v) 



VvCo 



b / ?i!(n')! 

2l'l+l''l(n+|/|)!(n'-f |/'|)! 

Jo Jo 



(63) 



2ii 



br' 



2\ (l'l + l''l)/2 



br^ 



= (^-tkovlB^^/2/b^ VyCWnl{n')l{n + |/|)!(n' + |/'|)! e-2(Mi,«)Vf' 
p!(-l)'=+'^'4'"'''(2fc§Z|i?V6) 



EE 

K=0 k'=0 



{n - K)l{n' - K')l{\l\ + Ky.{\l'\ + k'Y.kIk'I ' 



Here, we have defined 2p |/| + |/'| - |/ - l'\ + 2k + 2k', 
and made use of the series expansion for the Laguerre 
polynomial. 



The Hamiltonian matrix elements ( 62 1 and ( 63 ) can 



be solved numerically. We obtain a quantum dot energy 
spectrum for typical experimental conditions, as shown 
in Fig. |4j The device parameters are similar to those 
used in Fig. [T] Note that we have limited the analysis to 
orbitals in the range (2n -I- |/|) < 5. The different types 
of behavior observed in Fig. [4] are easily understood. In 
(c), the interface tilt is large, causing a strong suppres- 
sion of the valley coupling. Indeed, every energy level 
is doubly-degenerate. The apparent energy spectrum is 
equivalent to a circular dot in a single- valley material. In 
(b), a smaller, but physically realistic, interfacial tilt lifts 
the orbital and valley degeneracies, and causes a broad- 
ening of the high-field energy "bands." In this case, the 
tilt angle is still large enough that the valley splitting 
is smaller than the orbital splitting at zero field. We 
observe a complex assortment of level crossings and anti- 
crossings. In (a), the interfacial tilt is small enough that 



the valley splitting is larger than the orbital splitting at 
zero field. At high fields, the levels separate into two 
distinct subbands. 

The emergence of well-defined valley states is not ubiq- 
uitous in Fig. |4] By "well-defined" here, we mean that 
several states have identical valley characteristics, so that 
the valley index remains a good quantum number, even 
after the perturbation. Well-defined valley states are un- 
common because valley-orbit coupling mixes the unper- 
turbed orbitals differently for each eigenstate. However, 
when inter-orbital valley coupling is weak, well-defined 
valley characteristics can emerge. 

Weak inter-orbital valley coupling typically occurs in 
two ways. First, it can occur when the valley splitting 
is inherently small compared to the orbital splitting, for 
example when the interface tilt is large. In this case, the 
intra-orbital valley coupling (i.e., between states with the 
same orbital quantum numbers) is small, and the inter- 
orbital valley coupling is much smaller. The valley po- 
tential then produces weak pair-wise couplings between 
states with the same orbital quantum numbers, as in 
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2 4 6 8 

Magnetic field, 6 (T) 

FIG. 4: Energy levels for a circular quantum dot in a per- 
pendicular field geometry, including valley coupling. Dot 
parameters are appropriate for a single-electron spin qubit: 
hujo = 0.75 meV, v^Co = 0-5 meV. Uniform step orientation 
is set to ^ = 7r/4. The effective miscut tilt angles are given 
by (a) i9 = 0.03°, (b) ■& = 0.3°, and (c) i9 = 3°. (b) is prob- 
ably closest to experimental conditions. In (c), each curve is 
doubly degenerate. 



splittings are comparable in size, as in Figs. |4]ja) and 
(b), the inter-orbital coupling competes with the intra- 
orbital coupling, leading to a more complicated valley- 
orbit mixing. 

Well-defined valley states also emerge in the high-field 
limit, as observed in Fig. |4]^a). Within a given band, we 
find that valley coupling occurs primarily between states 
with the same angular momentum quantum number, and 
that the valley splitting is identical for all such pairs. For 
states with different I, Eq. ( [63| shows that the valley cou- 
pling is proportional to ', which vanishes at 
high fields. Thus, in the n = manifold, we observe pairs 
of states with the energy splitting 2A077, with many in- 
tervening states between them. In other words, two sub- 
bands emerge, with a uniform splitting. The ground state 
gap between the two lowest levels is then determined by 
orbital effects, rather than the valley coupling. Within a 
given subband, the eigenstates all have an identical valley 
character: one subband is "even," the other is "odd." 

In spite of the rather complicated behavior observed 
in Fig. |4][a)-(c), a simple expression can be given for the 
splitting between the two lowest energy levels. We first 
note that the two lowest states in Fig. |4|^c) are nearly 
degenerate, with orbital quantum numbers n = I — 0. In 
Fig. |4][b), the low- lying states both have an n = / = 
character at low fields. However, the higher state crosses 
over to an n = 0, / = — 1 character at high fields. In 
Fig. |4|^a), the first excited state has an 71 = 0, ^ = — 1 
character over the entire field range. From such consid- 
erations, we see that the ground state gap should be well 
characterized within the manifold spanned by quantum 
numbers n = 0, I = {—1, 0}, and v = ±1. Evaluating the 
perturbation Hamiltonian in this manifold again leads to 
Eq. (46), where we now use the definitions (61 1, (65), and 



(66). The results for the ground state gap, Eq. (49), and 



the mixing fraction, Eq. (50), then apply, if we make the 
replacement 



hiOy — > {hiOc/2){b— 1). 



(67) 



Fig. ^c). In this case, the valley splitting is different 
for each pair of states. For example, in the n — mani- 
fold (the lowest "band" in the high-field regime) , we find 
that 



(64) 



which depends on the quantum number I. Here, Ag is 
the same as before, but 



(65) 
(66) 



?7 = e ^ 



and Ln{x) = L^{x) is an ordinary Laguerre polynomial. 
The valley eigenstates have an identical character for all 
orbital states, however, and can be classified as either 
"even" or "odd." Similar behavior is observed for large 
interfacial tilt angles in Fig. [l] When valley and orbital 



We emphasize that these results correctly interpolate be- 
tween orbital and valley dominated behavior, over the 
entire field range. 



B. Elliptical Dot 

To close this section, we investigate an elliptical quan- 
tum dot in a perpendicular magnetic field. We cannot 
directly apply the results of Sec. |III[ However, the prob- 
lem is still separable with respect to the z' coordinate, 
as in Eq. (59). The magnetic field now couples the two 

and y', in a fundamental way. 



transverse directions, x' 
The transverse effective mass is nearly isotropic when the 
interfacial tilt angle d is small.'^ rux — ruy ^ mt. So it 
is unnecessary to rescale the length axes. Dropping the 
prime notation, the transverse envelope function is given 
by 
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iJxy MiLOix^ + 2iLMiujiM2UJ2xy + M2Uj2y'^ 



Cfc ; 



2h{L^MiUJiM20J2 + 1) 
y/2MiU!i/h{x + iLM2UJ2y) 



(68) 



fc=0 1=0 



k\l\{n^ - k)\{ny - l)\ 



L2M1W1M2W2 + 1 



yj2M2L02 / f^{LMiLOiX — iy) 
L'^MiUJiM20J2 + 1 



where, now, 



mt 

2uir 



(69) 



Wl9 



(72) 



±y^(w2+^2+^2)2_4^2^2 



-11/2 



L = 



"m-t^J {i^l + w2 + ^2)2 _ AuIujI 



(70) 



Mi,2 = (71) 
2mt^{ujl + a;2 + a;2)2 - 4a;2^2 



(a;2 - a;2 ± ^2) _,_ y^(cj2 + ^2 + ^2)2 _ 4^2^ 



rig; ^riy 



(-z)"« (3/4)1/4 



(na.)!(ny)!VMia;iM2a;2 
;i(L2Mia;iM2W2 + 1) 



1/2 



(73) 



(_l)-V2i2a+^-)/4(3 _ 4^2G2)0-.-i)/4(i _ 2i^2)(.-o/2pg_+;+i)/2 (^^|^) (fc + Z even), 
(fc + ^odd), 



G 



2L'^MiujiM20J2 
L2M1CJ1M2W2 + 1' 



L2Mia;iM2W2 + 1 ' 



(75) 



(76) 



Note that and riy are non- negative integers, P^fiz) is 
an associated Legendre function, and we have assumed 

The energy eigenvalues are given by 



1 



-■n^.ny — 2 1^2/ ^ \"'y • 2 



1 



As before, we have adopted the lowest subband approxi- 
mation. 

We can compute the valley coupling matrix elements. 
The intra-valley terms are given by 



{n'xn'yV\H\na,nyv) = 5„,,<(5„„,„^ [en.,n„ + VvC.1], (78) 



+ [nx + -] huji + [ ny + - \ huj2- (77) as in Eq. (10 1. The inter- valley terms are given by 
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3K)!K)!(n,)!(n;)! 



(79) 



X exp 



cos^ ip sin^ (fi 
AIiUJi M2UJ2 



«x < „,-fcn„-;n;-fc'n;-/' Lm/2J Lg/2J Lni72j L972J 

EEEEEEEEEEEE 

fc=0 (=0 k'=0 l'=0 m=0 (}=0 m'=0 i3'=0 p=0 r=0 p'=0 r'=0 



g + 1 



t/2 



xH, (koi) cos ip^/h{Q + 1)/Micji) i7t_, (^Msin(^y^i(QTl07M2^ 
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where we have defined 



spends to £0,0 hi Eq. ( 77 ) , and 



m + m' + q + q' — 2{p + p' + r + r') , 
Tlx + n'^ + riy + n'y — k — k' — I — I' 
-2{p + p' + r + r'), 
Tlx + n'x — k — k' — m — m' + q + q' 
-2(r + r'), 



and 



Q 



uj'^MiU}iM2UJ2/m^ 



(80) 
(81) 
(82) 

(83) 



and we have made use of the series expansion for Hermite 
polynomials!^ 

The valley coupling for an elliptical quantum dot in a 
perpendicular field geometry, Eq. ( 79 ) , is the most gen- 



eral solution obtained in this paper. The resulting energy 
spectra are qualitatively similar to the geometries studied 
above, and we do not plot them here. 

We consider two limiting cases. In the small dot limit, 
the ground state can be treated using the same pertur- 
bation Hamiltonian as Eq. ( 40 1 . However, Eq now corre- 



T] = exp 



hk^r{Q + i) 



cos^ f sin^ If 

MlLUl M2LO2 



(84) 



while Aq is unchanged. The ground state gap is then 
given by 



Eg ~ 2Aor/. 



(85) 



In the further limit of uj^ — ^y, we recover the equivalent 
result for the circular quantum dot in a perpendicular 
field geometry: 



Eg^2v,CleM~'^{kolB-&f/b\ 



(86) 



Alternatively, keeping uJx 7^ Wj,, but taking the limit B — > 
0, we recover Eq. (44 1, corresponding to a small elliptical 



quantum dot in the parallel field geometry. 

We also consider the important case where the valley 
coupling is of the same order as the lowest orbital split 
ting (Ao?? 



HID 



we 



hujy < fujjx)- Analogous to Sec. 
compute the Hamiltonian in the manifold spanned by 
the ordered basis {|00+), |00-), |01-|-), |01-)}, obtaining 



H ■ 



( £0,0 + Ao Aqt? 
Aof? £0,0 + Ao 

Ao?7(r + ig) 
VAo77(-r - iq) 





Ao?7(r - iq) 
£0,1 + Ao 



Ao77(-r + iq) 


Ao7?(l - r2 - g2) 



(87) 



Ao77(l-H-g^) £0,1 +Ao 



r 



where the various parameters are defined above, except 



(k^d sin (p) y/2h/M2UJ2, 
{ko'd cos ip) \/2hQlM\Ltj\. 



(89) 



We can diagonalize the Hamiltonian and compute the 
ground state gap and the mixing fraction as before. This 
yields identical results to Eqs. ( 49 1 and ( 50 1 , if we replace 
q^ — > ((?^ + T"^), and use the parameter definitions given 
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above. In the further hmit i? — )• 0, we find that r = 0, 
and the parameters defined in this section reduce to those 
of Sec. |IIIE)| The results are then equivalent. 



V. DISCUSSION AND CONCLUSIONS 

The term "valley splitting" is widely used in the litera- 
ture of indirect gap semiconductors. However, its mean- 
ing is not precise. As demonstrated above, valley-orbit 
coupling can lead to a complete mixing of orbital wave- 
functions in the regime of interest for quantum comput- 
ing, so that pure valley and orbital quantum number are 
no longer good quantum numbers. When atomic steps 
are present at the quantum well interface, many different 
"valley splittings" can be observed in the energy spec- 
trum, none of which are universal. 

We propose some formal definitions to help clarify this 
situation. First, we define the "characteristic" or the- 
oretical maximum of the valley splitting, corresponding 
to 2Ao. This splitting is achieved in a hypothetical flat 
quantum well, and its value is independent of the quan- 
tum dot shape. Second, we define a "renormalized" val- 
ley splitting, corresponding to 2Ao?7. This quantity ac- 
counts for the suppression caused by a tilted interface. 
The renormalization factor rj depends on the dot shape, 
but it does not include the effects of inter-orbital valley 
coupling. Finally, we note that the energy splitting be- 
tween the lowest two quantum dot eigenstates should not, 
accurately, be referred to as the valley splitting, since this 
usage is only valid in limiting cases. We have identified 
some of those cases here {e.g., the small dot limit and the 
high-field limit) . More accurately, we should simply refer 
to the "ground state gap," which is of particular interest 
for quantum computing and other applications. In this 
paper, we have obtained analytical expressions for the 
ground state gap and for the mixing fraction, in several 
cases of interest. 

Our main conclusions for the valley coupling in quan- 
tum dots can be summarized as follows: the shape and 
the size of the quantum dot, the number of steps that it 
covers at the quantum well interface, and the nature of 
the disorder all affect the suppression of the valley split- 
ting. The trends in the energy spectrum of Fig. [l] are dif- 
ficult to characterize in a simple way. One might expect 
excited orbitals to exhibit a larger suppression because 
of their larger size. However, the effect is obscured by 
the nodal structure of the wavef unctions. The behavior 
is more plain when we compare ground state wavefunc- 
tions of different sizes, as in Fig. [3] In this case, smaller 
dots clearly exhibit a larger valley coupling. As expected, 
the effect depends sensitively on the tilt angle of the in- 
terface. The dependence of the valley coupling on dot 
size is also manifested in Fig. |4|^a), at high fields. Here, 
the interface tilt angle is rather small, and the magnetic 
confinement of the quantum dot is strong, so tilting has 
little effect on the energy spectrum. The valley splitting 
then approaches its theoretical maximum, regardless of 



the specific dot shape. 

Experiments in Si 2DEGs generally corroborate our re- 
sults. Examples include Hall bar measurements, where 
the valley splitting grows monotonically with the mag- 
netic field. ■^^ This behavior is consistent with the 
case of small or vanishing electrostatic confinement in a 
quantum dot [e.g., the case Z = 0, & ~ 1 in Eq. (64|]. In 
a quantum point contact, valley splitting is significantly 
enhanced due to lateral confinement.'^ In a Si/SiGe two- 
electron quantum dot, a large singlet-triplet splitting be- 
tween 0.1 and 0.3 meV has been measured at zero-field.!^ 
Since the triplet state involves an excited single-electron 
orbital, the singlet-triplet splitting provides a good esti- 
mate of the ground state gap. This solitary data point is 
therefore consistent with Fig. |4|^b). Due to valley-orbit 
coupling, it is not appropriate to identify the splitting as 
purely orbital-like or valley-like. However, small gaps of 
any type are anathema for qubit operations. The large 
value of the measured gap implies an acceptable splitting 
for qubit applications.^"* In future experiments, it would 
be desirable to perform excited state spectroscopy in sim- 
ilar dots, as a function of magnetic field. The present 
theory would then allow the valley coupling and other 
quantum dot confinement parameters to be determined. 

In this paper, we have paid special attention to the 
inter-orbital valley coupling. The origin of this coupling 
can be understood in terms of the dot shape. According 
to the variational principle, a quantum dot will minimize 
its ground state energy by maximizing its valley split- 
ting. In our perturbation theory, the ground state will 
therefore mix with asymmetric excited states that allow 
it to squeeze along the direction of the step gradient. The 
mixing occurs only between opposing valley states, v and 
V, whose coupling is the source of valley splitting. 

The inter-orbital component of valley coupling has not 
been emphasized previously. In quantum information 
applications, it leads to a potential source of decoher- 
ence for spin qubits. When phonons are present, valley 
scattering can occur in a process analogous to spin flip 
transitions!^ Two separ ate ingredients are needed to flip 
a spin by a phononl^^'SS spin-orbit coupling (to mix spin 
states via the excited orbitals) and a magnetic field (to 
lift the Kramers degeneracy) . Valley scattering can occur 
in the presence of two analogous processes: valley-orbit 
coupling at a tilted interface (to mix valley states via the 
excited orbitals) and a sharp confinement potential (to 
lift the valley degeneracy). 

In the context of quantum computing, valley scatter- 
ing does not directly affect the spin state or lead to spin 
decoherence. However, there is an indirect spin effect, 
mediated by the Pauli principle, which depends on the 
valley coupling. The predominant two-spin gate inter- 
action is the exchange coupling,^ which occurs because 
any two electrons must have different quantum numbers. 
When the valley degree of freedom is introduced, the size 
of the Hilbert space doubles, providing an opportunity 
for an electron to "leak" into a non-qubit sector via val- 
ley scattering. If this happens, the electron will have a 
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different quantum number, and the exchange couphng 
will be ineffective. The spin state is then uncontrollable. 
Aside from this mechanism, a direct spin-valley coupling 
has also been predicted for Si 2DEGs, with a magnitude 
comparable to the spin-orbit coupling.l^^'^ We do not 
consider this effect here. 

Finally, we note that the results obtained in this paper 
were obtained by assuming a smooth, uniformly tilted 
interface. We expect such results to pertain to more gen- 
eral situations. However, step disorder and discreteness 
can both have a quantitative effect on the valley split- 
ting. In the Appendix, we discuss several experimental 
conditions that are inconsistent with our approximations, 
including: (i) step discreteness, which becomes impor- 
tant for widely separated steps; inhomogeneous tilting, 
which affects (ii) the quantum dot confinement potential 
and (iii) the delicate phenomenon of destructive interfer- 
ence in the valley coupling; (iv) step bunching, which also 
tends to moderate the destructive interference. Since all 
of these phenomena enhance the valley splitting, we can 
treat them, phenomenologically, through an effective tilt 
angle that is smaller than the nominal tilt. 

In the absence of any detailed knowledge of the interfa- 
cial disorder, 'd must remain a phenomenological param- 
eter. However, the relation between the average tilt and 
the effective tilt, and the dependence on different forms 
of disorder, remains an important topic for future inves- 
tigation. For example, the energy levels of a Si double 
quantum dot could potentially be measured as a function 
of magnetic field by the method of Ref. Compari- 
son with the present theory can then provide information 
about More intriguingly perhaps, it may be possible 
to use gating methods to shift the center of a dot, and 
thereby map out as a function of position. 

In conclusion, we have developed a theory of valley 
coupling in realistic quantum dots. The resulting energy 
spectra exhibit crossings and anti-crossings, as a function 
of the interfacial tilt angle and magnetic field. Due to 
valley-orbit coupling, the ground state gap is not strictly 
orbital-like or valley-like, except in certain limiting cases. 
For quantum dots of interest in quantum computing, we 
find that inter-valley orbital coupling plays an important 
role in device operation, and in determining the ground 
state gap. 
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Appendix A: Validity of Approximations Used in 
the Calculations 

This appendix addresses the approximations made in 
the main text that enable analytic solutions to be ob- 
tained, and it discusses the conditions under which they 
apply. The main approximations are that: (1) the inter- 
facial tilt can be treated as continuous and the presence 
of individual atomic steps can be ignored, and (2) the 
disorder in the step locations can be ignored. The first 
approximation suggests that there will be no relaxation 
of the electron wavefunction to conform to the pattern 
of steps at the interface. The second approximation is 
equivalent to assuming the interfacial tilt is smooth and 
uniform. 



1. Step Discreteness: Phase Relaxation at Discrete 

Steps 

In this work, we treat the interface as continuous 
and we ignore the presence of individual atomic steps. 
Naively, one would expect such an approximation to be 
reasonable because the step is much smaller than the 
characteristic size of the envelope function. For exam- 
ple, the height of a single atomic step along [001] is 
a/4 = 1.4 A, while a quantum dot may extend over lO's 
of nm. 

Energy arguments suggest that such an argument 
breaks down when the interface tilt angle is small. In 
Eq. ^, we have taken the complex phase difference be- 
tween an^n,v a-nd an,n,v to be a constant, independent of 
the lateral position. The global value of the phase differ- 
ence is determined by minimizing the ground state energy 
for the whole quantum dot, according to degenerate per- 
turbation theory. Since the optimal value of the phase 
difference depends on the interface position zq , it is clear 
that valley coupling cannot be optimized simultaneously 
on both sides of a step while the phase remains constant. 
By allowing the phase to vary, we can reduce this valley 
coupling energy penalty, but we pay a price in kinetic 
energy. Careful energy balance calculations show that 
phase relaxation becomes favorable when -d < 0.1 — 0.3°. 
In the present work, phase relaxation would tend to re- 
duce valley-orbit mixing. 



2. Disorder-Induced Confinement Potentials 

Step disorder causes non-uniform tilting of the inter- 
face, and leads to spatial variations of the valley splitting 
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and the quantum dot confinement potential. The posi- 
tion of the electron will shift when the local curvature 
of the valley splitting is larger than the curvature of the 
electrostatic potential. The net effect of the shift is to 
reduce the effective tilt angle. 



single-atom steps are essentially a worst-case scenario, 
leading to maximal suppression of the valley splitting. 
The dependence of the effective tilt angle on step bunch- 
ing is therefore rather complicated. 



3. Disorder-Induced Symmetry Reduction 



5. Summary 



For highly uniform step distributions, the suppression 
of the valley coupling by valley phase interference is very 
effective. However, this interference effect is somewhat 
delicate. The reduction of symmetry by step disorder 
can lead to an order of magnitude increase in the val- 
ley splitting.l^S Related effects are also observed in the 
presence of alloy disorder in the SiGe quantum well 
barriers.l^ While step disorder tends to reduce the effec- 
tive tilt angle,!^ the net effect of alloy disorder appears 
to be more complicatedi ^^ ' ^^ ' 



4. Step Bunching 

Interfacial steps may become bunched due to random- 
ness, or as a consequence of strained growth.'SSl The sup- 
pression of valley splitting, as discussed in this paper, is 
particularly effective for single-atom steps, as compared 
to bunched steps. We can understand this as follows. 
From Eq. ^ and the related discussion, the phase dif- 
ference between the z valleys for a flat interface is 2fco2o- 
At a step of height a/4, we see that zq ZQ + a/4, so the 
phase difference suddenly changes by 0.827r. Since the 
two sides of the step are almost fully out of phase, there 
is a significant suppression of the total valley splitting. 
On the other hand, a two-fold bunched step of height 
a/2 has a much milder effect on the valley coupling. In 
fact, the first step geometry with a stronger effect on the 
valley coupling is a six-fold bunched step. In this sense. 



To summarize, the assumption that the steps can be 
treated as continuous is expected to be valid, but val- 
ley phase relaxation and variability in the step density 
both serve to modify the effective tilt angle of the inter- 
face. Assuming that an appropriate effective tilt angle 
is used, we expect our analytical expressions to provide 
an accurate depiction of the valley coupling for typical 
quantum dot geometries. Our results therefore provide 
a useful first approximation for understanding the valley 
coup ling. Se veral numerical results confirm this point of 
view.™llll 

To obtain a quantitative estimate for the valley cou- 
pling, one must know the effective interfacial tilt angle "d. 
For the reasons described above, the effective tilt depends 
in a nontrivial way on the average tilt angle, the disorder 
in the step configuration, and the materials properties 
of the interface. A quantitative calculation of is par- 
ticularly challenging, since the dependence of the valley 
coupling on the tilt angle is exponential, so that errors 
and fiuctuations in the tilt angle are essentially magni- 
fied. For example, in Eq. (451, an exponential suppres- 



sion of the valley coupling was predicted for tilt angles 
larger than 0.26°, based on realistic dot parameters. Ex- 
perimentally, a large valley splitting has been observed 
for wavefunctions in a quantum point contact with a 2° 
misculpi' and in quantum dots in samples with nominal 
tilt angles approaching 2° (Ref. - evidence that the 
effective tilt in these devices is significantly smaller than 
the nominal tilt. 
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